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Abstract 

o 

. A path integral measure for gravity should also preserve the fundamental symmetry of 

general relativity, which is diffeomorphism symmetry. In previous work, we argued that 
a successful implementation of this symmetry into discrete quantum gravity models would 
imply discretization independence. We therefore consider the requirement of triangulation 
independence for the measure in (linearized) Regge calculus, which is a discrete model for 
I quantum gravity, appearing in the semi-classical limit of spin foam models. To this end we 

develop a technique to evaluate the linearized Regge action associated to Pachner moves in 
3D and 4D and show that it has a simple, factorized structure. We succeed in finding a 
^ . local measure for 3D (linearized) Regge calculus that leads to triangulation independence. 

I ' This measure factor coincides with the asymptotics of the Ponzano Regge Model, a 3D spin 

. foam model for gravity. We furthermore discuss to which extent one can find a triangulation 

independent measure for 4D Regge calculus and how such a measure would be related to a 
quantum model for 4D flat space. To this end, we also determine the dependence of classical 
Regge calculus on the choice of triangulation in 3D and 4D. 

^ ; 1 Introduction 

\o 

Many approaches to quantum gravity, such as spin foams [T], group field theories [2], (causal) 
dynamical triangulations [31 0] and Regge quantum gravity [5], rely on a path integral approach. 
A (non-perturbative) path integral has to be regularized to make it well defined. In the process 
of this regularization, several choices have to be made, that differ in the various approaches. 
Broadly one can understand these choices as deciding on a measure on the space of all geome- 
^ ■ tries. This includes various aspects, such as to define the space of geometries, for example the 

space of all triangulations with fixed edge lengths in dynamical triangulations versus the space 
defined by allowing all possible edge lengths (satisfying generalized triangle inequalities) in a 
fixed triangulation such as in Regge calculus, or some generalized discrete geometric spaces, as 
appearing in loop quantum gravity [U [7] . A related question is whether to include a sum over 
triangulations, such as in (causal) dynamical triangulations and group field theories, or even 
over two-complexes as suggested for spin foam models [5]. Alternatively the path integral may 
just include an integration over geometric data associated to a given, fixed, discretization. For 
discussions on the relation between these approaches, see ITU| [TT]. 

One reason for this many different suggestions, is that the space of all (discrete) geometries 
and its relation to the corresponding continuum space needs to be better understood [IH [13] . 
Many difficulties are rooted in the role of diffeomorphism symmetry, by which the space of 
metrics has to be quotiened to obtain the space of geometries. Discretizations obscure the 
role of diffeomorphisms, see |10^ [H] for a discussion. In particular, for a precise notion of 
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diffeomorphism symmetry in the discrete [10], one can show that this symmetry is broken for 
4D Regge gravity [15]. However, if this symmetry would hold in discrete gravity, one could 
hope for a unique anomaly-free (with repsect to diffeomorphisms) measure [16]. As is also 
argued in |16l I17j . the implementation of this symmetry into discrete gravity (of Regge type, 
i.e. with geometric data on a fixed triangulation or discretization), would make such a theory 
triangulation or discretization independent. In this case there would also be no need of summing 
over triangulations, which is often employed to obtain a triangulation independent theory. 

One can expect to find such a discretization independent theory for 3D gravity, which is a 
topological theory, i.e. there are no local physical (propagating) degrees of freedom. In fact 
we will succeed to find a triangulation invariant path integral description for 3D (linearized) 
Regge calculus. 4D gravity features local propagating degrees of freedom and a discretization 
independent model will require a non-local structure and moreover control over the solutions of 
the system [18]. Nevertheless, as argued in |16j the choice of path integral measure is important 
for the convergence of the model, also under a renormalization flow, which might be employed 
to find improved discretizations [19]. Moreover 4D classical Regge gravity is invariant under a 
set of certain local changes of the triangulation. One might therefore ask also for invariance of 
the path integral under this set of local changes. 

In this work we will concentrate on finding a measure in a (Euclidean) Regge calculus set up, 
that is as much triangulation independent, as possible. Before explaining this in more detail we 
will shortly review different measures suggested so far in the literature |20j . One method would 
be to discretize the (formal) continuum path integral 

Icont= Yi dgpr{x) ]^(^Jdet{g^)j exp{-SEH) ■ (1-1) 

X,p>T X ^ ^ 

Here Seh is the (Euclidean) Einstein Hilbert continuum action and (y^det((7^i/))" is a factor 
which can be obtained from the DeWitt metric on (geometric) superspace [21]. More specifically 
the DeWitt measure [21] prescribes a = in 4D and a = — 1 in 3D. However also other values 
of a have been suggested [20], for instance a = — (D + 1) for the Misner measure [22], where D 
is the dimension of space time. A priori it is not clear which choice to prefer |20j . 

Regge calculus [5] provides a discretization of the Einstein Hilbert action Seh-, defined 
on a triangulation. The metric data are replaced by edge lengths 1^ associated to the edges of 
the triangulation. As \J det((7^^) gives the local space time volume a natural discretization of 
this factor is given by the volumes Va of the top-dimensional simplices A, i.e. 4-simplices in 
4D and tetrahedra in 3D. A straightforward discretization of (jl.ip is then given by [20] (modulo 
numerical constants) 

W= /H'^^'n^A exp(-S^j) . (1.2) 

Concerning the range of integration it will always be understood that the generalized triangle 
inequalities are satisfied. These require positive volume for all (sub-) simplices and are therefore 
equivalent to restricting the integration in the continuum path integral (jl.ip to positive definite 
metrics. Apart from this requirement of triangle inequalities (which are technically very difficult 
to implement) the measure used in (|1.2|) has the advantage of being especially simple, in particular 
localjj The simplicity is also a reason why a = in 4D seems to be preferred |20] . 

^Another suggestion is to use a measure of the form W^V^^ dl^, which is scale invariant. (The Regge action 
without cosmological constant term is invariant -up to an overall factor- under global rescaling of the edge lengths.) 
However, this measure did not lead to satisfying results in numerical simulations, see [13] and references therein. 
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In this paper we will consider a requirement of triangulation independence for the path inte- 
gral measure. This requirement is also connected [16] with a discrete notion of diffeomorphism 
invariance [T7]. Hence asking for triangulation independence amounts to requiring an anomaly 
free measure with respect to the diffeomorphisms, see also [23] for a discussion in the spin foam 
context. 

Specifically we ask for invariance of the (linearized) model defined by (|1.2p under Pachner 
moves [23]. These are local changes of the triangulation, that act ergodically, i. e. two topo- 
logically equivalent triangulations can always be transformed into each other by a sequence of 
Pachner moves. Restricting the measure to the local ansatz (|1.2|) we will find that our results 
suggest to fix the parameter a to a = — ^ both in 3D and in 4D. Interestingly this conforms 
completely with the semi-classical analysis \25\ [26] of the Ponzano-Regge model [27] in 3D. 
This is a triangulation independent (spin foam) model for 3D quantum gravity, based on dis- 
crete variables. The case of 4D is much more involved. Firstly, being an interacting theory 
with propagating degrees of freedom, one cannot expect to obtain a triangulation independent 
model, with just local interactions, as in the Regge action [18]. Indeed, we will precisely show 
in which sense the (linearized) 4D Regge action fails to be triangulation independent. Although 
the semiclassical analysis of the 4D models [281 [29] could show that the Regge action appears in 
a /i — )• limit of the amplitudes, the corresponding measure factor has not been specified yet as 
a function of the geometric variables. For future work it will be interesting to compare in more 
detail the spin foam results with Regge gravity. Also a measure ambiguity shows up in choosing 
so called edge and face amplitudes |23l [8l [30] . These ambiguities could also be restricted by 
asking for as much triangulation independence as possible, similar to the method proposed here. 
Hence it would be very interesting to study the behavior of spin foam amplitudes under Pachner 
moves [3T] . 

There are also other suggestions for the Regge measure, which are non-local. As these are 
far more complicated explicit computations they have mostly been restricted to 2D. The Regge- 
Lund measure |32^ [331 [20] is obtained by discretizing first the deWitt super metric and then 
taking the determinant (whereas in ()1.2p this is performed the other way around). The result is 
given by 



where 




and /3 is another ambiguous parameter. As the determinant has to be taken of a matrix, which 
is indexed by all the edges of the triangulation, the result is potentially quite non-local. Further 
discussion of this measure can be found in [331 [2D] . 

In 3D, where gravity is a topological theory, we will find that a local measure is sufficient 
to guarantee triangulation independence of the (linearized) theory. In 4D, as previously men- 
tioned one cannot expect to find complete triangulation independence for the path integral as 
already the action is not triangulation independent. (More precisely it is the Hamilton- J acobi 
function as a functional of the boundary data, that is not invariant under the change of the bulk 
triangulation.) One can however ask for invariance under a restricted set of Pachner moves, 
under which the action happens to be invariant. These are the 4 — 2 and 5 — 1 moves (but 
not the 3 — 3 move). Nevertheless also for these moves we will find that a factor appears that 
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features a certain non-local structure. At this stage it seems however more promising to con- 
struct improved measures and actions directly by coarse graining and the method of perfect 
discretizations [19 } [THl [T6] . 

Ultimately, another criterion that any quantum gravity model has to satisfy, is to display 
the correct large scale limit. Also here a measure term could be essential. For investigations in 
2D Regge see for instance |34j . for discussion of the influence of the measure in the context of 
dynamical triangulations see [35\ IBB] . Another suggestion for constructing a measure for Regge 
gravity, is to mod out a certain subgroup of the continuum diffeomorphism group |37j . This 
results again in a highly non-local measure, where explicit results are mostly restricted to 2D. 

In the next section we introduce the Regge action and its expansion up to second order. 
This requires the calculation of its Hessian matrix, which will be one of the main subjects of this 
work. Furthermore we discuss the concept of Pachner moves and briefly present the Pachner 
moves in 3D and 4D. Section [3] deals with a general method to compute the Hessian matrix in 
3D and presents the application of this method to the Regge actions associated to the Pachner 
moves. Then we examine invariance of the path integral under Pachner moves and define a 
suitable measure factor. The results for 3D will be summarized in section 21 In section [5] we 
extend our method to compute the Hessian matrix in 4D, examine invariance of the path integral 
under Pachner moves and discuss a suitable measure. The results in 4D are then summarized 
in section [BJ We conclude this work with a discussion of our results in section [71 



2 Linearized Regge Calculus 

The Regge action (which we will denote by S in the following) provides a discretization of the 
Einstein Hilbert action for gravity. It is defined on a triangulation, the geometry is a piecewise 
flat one, and the geometric data are encoded in the lengths of the edges in this triangulation. 
For Regge type actions based on different geometric variables nearer to spin foams, specifically 
areas and angles, see [38| [39]. 

In the following we will consider the Euclidean path integral for the Regge discretization of 
gravity on a given 3D or 4D triangulation 

f H dleKk)exp{-S} . (2.1) 

Jle\ecaM eCbulk 

Here le\e(idM denotes the boundary conditions, which we take to be fixed length variables for 
the edges in the boundary triangulation. /u(/e) is a suitable measure factor. In (j2.ip not all 
edge lengths combinations are allowed since the edge lengths have to satisfy generalized triangle 
inequalities, i.e. all the (2D, 3D and, in case, 4D) volumes have to be positive. The (Euclidean) 
Regge action in arbitrary dimension D can be written in the following form 

/iCbulk \ (jOd/i / hchdrj \ a^^h / 

where denotes D-simplices, i.e. D-dimensional simplices with D + 1 vertices, h denotes 
'hinges', i.e. D — 2-simplices, Vh is the volume of a hinge and 9^ ^ denotes the internal dihedral 
angle in the D-simplex at the hinge h. The terms in brackets in ()2.2p define the bulk and 
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boundary deficit angles 



(^"^■^^ := 27r - (2.3) 



where k depends on the number of pieces one is glueing together at this boundary. If there are 
only two pieces we have k = 1. 

The dihedral angles are complicated functions of the lengths variables, so that the integral 
in ()2.ip cannot be computed analytically. Additionally one has to take the generalized triangle 
inequalities for the range of integration into account. 

To circumvent this issue, we consider linearized Regge Calculus in which one chooses a 
classical background solution (for the edge lengths) if'^ satisfying the triangle inequalities and 
one quantizes, i.e. integrates over, the perturbations Ag around it. 

Therefore consider a small perturbation around a background solution 

le = if^ + Xe (2.5) 

and expand the Regge action up to second order in the perturbation variables Xe'- 

dS 



S = 



/g — jSi d\i f> 



1 d's 

iM"^ ^ ^ 2 dLdU 



I'e — i-e 



XeXe' ■ (2.6) 



The background edge lengths l'^^ are defined as the solution to the Regge equations, 

dS ^dVh 

^ toe ^ 

such that the first order term in (|2.6|) vanishes for the bulk edges. More specifically we take 
the background solution to be (locally) flat, that is o;^^"''^^ = 0. (This is exactly the equation of 
motion in 3D.) The second order term is defined by the matrix of second derivatives, that is the 
Hessian. 

In three dimension one obtains due to the Schliifli identit}!! 



(21 



dledle' die' 

In four dimensions we obtain (using again the Schlafli identity): 

d'^S sr^ dAhduh _ sr^ d'^Ah (bulk) _ d'^Ah (bdry) 

2^ Pll, m.. 2^ rlL.BL,^f^ 2^ rlLrlL,^'^ ' ^ ' 



dhdlpi ^ dip' dip ^ dlpdlp' ^ OLdL' 

^ h ^ /iCbulk ^ hcbdry ^ 

The bulk deficit angles w^^"^'^^ vanish on a flat (background) solution. 

For the evaluation of these Hessian matrices, we will need the first derivatives of the dihedral 
angles with respect to the length variables. A formula valid for simplices of arbitrary dimension 



The Schliifii identity EhcA(«) = ensures that terms with second derivatives of the dihedral 



he 

angles vanish. 
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D can be found in |40j : 



+ cos{eki){cos{ekh)cos{6km) + cos{0ih)cos{6im))^ ■ (2.10) 

In (j'i.lOp 9ki denotes the dihedral angle (in a D-simplex) between the two D — 1-simplices formed 
without the vertices k and / respectively. Ihm is the length of the edge between vertices h and 
m. Vfi denotes the volume of the {D — 1) simplex formed without vertex h in the D-simplex. V 
is the volume of the respective D-simplex. 

In case the dihedral angle O^i and the edge Ihm do not share a vertex, which implies that 
(kl) = (hm), i.e. the hinge is formed without the vertices h and m in the D-simplex, equation 
(|2.10p simplifies using the convention cos On = — 1: 

dlhm D^sm Ohm y ^ 

1 hmyhym . p: 1 hmyhm 1 1 ^ 

■sm0/,™ = -— — — - — . (2.11) 



]j2 y2 D{D-l) V 

This result (j2.11|) will be crucial for an alternative derivation of the matrix elements of the 
Hessian ()2.1ip in section [3l This alternative derivation applies to configurations defining Pachner 
moves, which we will discuss in the next section. 



2.1 Pachner Moves 

Pachner moves are local changes of the triangulation which, if applied consecutively, allow to 
go from any triangulation of a given manifold to any other triangulation of that manifold [M] . 
In quantum Regge calculus one usually fixes the triangulation and just integrates over the 
edge lengths in this given triangulation. Given this definition the question arises of how the 
result depends on the choice of triangulation. Note that the triangulation is only an auxiliary 
structure, which is put in in order to regularize the (continuum) path integral. Hence it would be 
advantageous, if the path integral (with or without given boundary triangulation and condition) 
would depend minimally on the choice of (bulk) triangulation. In case the path integral does 
not depend at all on the triangulation, we do not even need to take any refinement limit (here 
of the bulk triangulation only), as the result will not change under refinement. Such a strong 
version of discretization independence can actually be expected in 3D, in which gravity is a 
topological theory, describing the dynamics of only global (topological) variables. Indeed we 
will find a measure that will render the path integral discretization independent in this sense. 
(That the linearized action is invariant under refinements has been shown in [IB].) Locally 
this discretization independence implies that the path integral is form invariant 'under Pachner 
moves'. More precisely, we will consider here Pachner moves arising by integrating out certain 
edges in the triangulation, so that the remaining edges still define a triangulation. The form of 
the (discretized) path integral should then be invariant. 

We will consider a similar requirement in 4D. This defines however, a theory with local 
degrees of freedom, where not even the (linearized) action is invariant under change of trian- 
gulation |4H [THl B2] . This broken invariance can however be isolated into one of the Pachner 
moves, the 3 — 3 move. Hence we can at least ask whether it is possible to define a measure that 
would render the path integral invariant under the remaining Pachner moves. 
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A X — y Pachner move changes a complex of x D-simplices into one with y L)-simphces 
without changing the boundary triangulation. Here the parameters x, y are related hjx + y = 
D + 2. Since the boundary is not changed, the Pachner moves act locally in the triangulation. 
This also allows us (for the cases with x > y) to consider the initial configuration oi x D 
simplices, to integrate out the bulk edges and to re-interpret the resulting partition function as 
one for the complex with y simplices. Note also that the x — y and the y — x move are inverse 
to each other. 

In the following sections we will introduce the Pachner moves in 3D and 4D and shortly 
point out some points pertaining to the dynamics defined by Regge calculus. 

2.2 Pachner moves in 3D 

Here we have two Pachner moves 3 — 2 and 4—1 (and their inverses) . Note that the equation 
of motion for 3D Regge calculus require flatness, i.e. vanishing deficit angles. 

2.2.1 3-2 move 

The first Pachner move we will consider is the 3 — 2 move, see Fig. [TJ In the initial config- 
uration three tetrahedra (0123), (0124) and (0134) share an edge (01). This is the only bulk 
edge. Removing (i.e. integrating out) this edge and introducing a triangle (123) we obtain a 
configuration of two tetrahedra (0234) and (1234) sharing this triangle. 

As there is only one bulk edge in the initial configuration, we will also have only one equation 
of motion. This equation of motion requires the vanishing of the bulk deficit angle ujqi and in 
this way fixes the length /qi of the edge (01) as a function of the boundary edge lengths. 








1 1 



Figure 1:3 — 2 move. The two tetrahedra can be split into three by connecting the two vertices 
separated by the shared triangle. The dashed line in the three tetrahedra configuration is the 
dynamical edge. 

2.2.2 4-1 move 

The other Pachner move in 3D we are going to discuss is the 4 — 1 move. Here in the initial 
configuration four tetrahedra share one vertex. This configuration can be obtained by subdi- 
viding a tetrahedron (1234) into four tetrahedra by placing one vertex into the tetrahedron 
and connecting with the other four vertices. In the 4 — 1 move this vertex and the adjacent 
edges are removed, leaving us with one tetrahedron (1234), see also Fig. [2j 
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In the initial configuration with four tetrahedra, there are four bulk edges, and hence four 
equations of motion. These, again require that the (four) bulk deficit angles have to vanish, i.e. 
that the complex has to be flat. We know that we can easily construct such solutions by placing 
a vertex into the (flat) tetrahedron (1234) and determining the lengths of the four bulk edges. 
There is, of course, a three-dimensional parameter space of where to place the inner vertex 
exactly, hence the solutions are not uniquely determined. This is the well known gauge freedom 
in Regge calculus on flat solutions [43 [ HOj [TU] , a discrete remnant of the difFeomorphism 

symmetry in the continuum. From this it follows that of the four equations of motions only one 
is independent and that we have to expect three null modes in the Hessian matrix of the system, 
signifying three gauge degrees of freedom. Further discussions and extensions to the case with 
cosmological constant can be found in [lOl [151 US]- 



4 4 



1 




2 2 



Figure 2: 4 — 1 move. The tetrahedron is split into four by placing one additional vertex inside 
the tetrahedron and connecting it to the remaining vertices in the boundary giving four internal 
edges (dashed). 



2.3 Pachner moves in 4D 

2.3.1 4-2 move 

This Pachner move is very similar to the 3 — 2 move in 3D, see Fig. [3l The initial configuration 
is one with four 4-simplices, (01234), (01235), (01245) and (01345), sharing one edge (01). All 
the other edges are boundary edges. Removing this edge and introducing a tetrahedron (2345), 
we obtain a configuration with two 4-simplices (02345) and (12345) sharing this tetrahedron. 

As there is only one bulk edge we again have only one equation of motion for the initial 
configuration. A (flat) solution can always be constructed in the following way: The boundary 
triangulation is the same as for two 4-simplices sharing one tetrahedron. Such a configuration 
can always (i.e. for all boundary edge lengths satisfying the appropriate inequalities) be em- 
bedded into flat 4D space. We can hence straightforwardly determine the distance between the 
vertices and 1 in the induced metric, which defines the length of the edge (01). 

In some exceptional cases there might be also solutions with curvature [15], however this 
seems to be rather a discretization artifact. For the perturbative solutions around flat space 
we are interested in, we can note that the linearized equations of motion have a unique (flat) 
solution for all (linearized) boundary perturbations. 

2.3.2 5-1 move 

The 5 — 1 is again analogous to the 4—1 move in 3D. In the initial configuration five 4-simplices 
share one vertex (0) which is adjacent to five bulk edges. Removing this vertex and the adjacent 
edges we are left with just one simplex (12345), see also Fig. [3| 
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Also here, we can construct for all boundary configurations flat solutions to the equations 
of motion. These can be found by placing the vertex (0) into the (flat) 4-simplex (12345) and 
determining the induced lengths of the edges (Ox), where x = 1, . . . , 5. For given boundary data 
there is a four-parameter space of such solutions, according to the four parameters describing 
the position of the vertex inside the 4-simplex. Hence we can expect four null modes for the 
Hessian of this configuration. 



1 1 




3 4 3 4 



Figure 4: 5 — 1 move. The 4-simplex is split into five 4-simplices by placing one vertex inside the 
4-simplex and connecting it to the boundary vertices, hence obtaining five bulk edges (dashed 
lines) . 



2.3.3 3-3 move 

We are left with the 3 — 3 move, which is significantly different from all the other Pachner moves 
discussed so far. 

Assume three 4-simplices (01234), (01235) and (01245) sharing one triangle (012). Note that 
this configuration does not include a triangle (345), as neither of the three 4-simplices contains 
the three vertices (3), (4), (5). 

The 3—3 move rebuilds this configuration into three 4-simplices (01345), (02345) and (12345) 
which share the triangle (345) and do not include the triangle (012), see also Fig. [5j 

In contrast to all other Pachner moves discussed so far the 3 — 3 move does not involve 
dynamical edges, i.e. all edges are in the boundary and therefore included in both configurations. 
We therefore do not have an equation of motion. Note however that, again in contrast to the 
other Pachner moves, not all boundary configurations define a flat geometry. That is, in both 
configurations we have only one bulk triangle. The vanishing of the deficit angle for this bulk 
triangle gives one condition for the length of the boundary edges. In case this condition is 
violated we do have a curved configuration. In particular, even on a flat background, we can 
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5 5 




14 14 



Figure 5:3 — 3 move. Three 4-simplices sharing the triangle (012) and not containing (345) are 
rebuilt into three 4-simpnces sharing the triangle (345) and not including triangle (012). The 
shared triangles are drawn dashed in this figure. Note that all edges are boundary edges and are 
contained in both configurations, so the configurations are determined by the shared triangle. 



have a curvature excitation, if the boundary perturbations do not satisfy the linearized flatness 
condition. 

In the following section we will specify the Hessian matrix of the Regge action associated 
to the various configurations appearing in the Pachner moves. We could start with the formula 
()2.10p for the derivatives of the dihedral angles to obtain the derivatives of the deficit angles, so 
that these can be combined to give the entries in the Hessian. This procedure would however 
result in very lengthy formulas and not use the fiatness of the background solution. We will use 
an alternative strategy, which will produce a quite enlightening structure for the Hessian, and 
for which we present some auxiliary formulas in the next section. 

3 Computation of the Hessian matrix in 3D 

In this section we will compute the matrix elements of the Hessian matrix. To do so one has to 
compute terms of the form ^ for which we will present a general strategy, similar to [461 ITT] . We 
need to extend the ideas in [46 ^ [TTj in order to also obtain the matrix elements of the Hessian 
indexed by edges in the boundary. First we will derive two auxiliary formulas, which is the 
subject of the next section. 

3.1 Auxiliary formulas 

For concreteness we will derive the auxiliary formulas for the initial configuration of the 3 — 2 
move (see Fig. [1]) with the bulk edge (01) as described in the previous section. Assume that 
this configuration can be embedded into fiat (3-dimensional) space, i.e. M^. This implies that 
for instance Iqi, the edge length of the dynamical edge (01), is fixed as a function of all other 
edge lengths. Hence there is one relation which all edge lenghts have to satisfy, which is wqi = 0, 
i.e. the deficit angle at the edge (01) vanishes. (Note that this relation can also be derived by 
requiring that the Cayley-Menger determinant of this configuration, giving the square of the 
4D volume, vanishes.) As there is one condition, at least two edge lenghts have to be varied in 
order to preserve this relation. Therefore consider variations of exactly two edge lenghts, I and 
I'. Alternatively one can interpret I' as a function of /, where all other edge lengths are (fixed) 
parameters. We will derive a relation between the edge variations 61 and 51' or alternatively the 
partial derivative 

For example, consider I' = I23 and / = l^^. We vary I23 and 1^4 such that wqi = 0, i.e. the 



10 



triangulation is still embeddable in flat space. This implies: 



= 5uj. 



01 



g{0123) 

^1 



59, 



(0124) 
01 



56, 



(0134) 
01 



where 0^^^^^^ is the dihedral angle at edge (01) in the tetrahedron (Olxy). 
neither depends on Z23 nor on ^34, as these edges are not part of the tetrahedron (0124). 
Using equation (|2.11|) . we obtain: 



.(0124) 
^1 



(3.1) 



since it 



89, 



(0123) 
01 



dl 



23 



^01^23 
6^1 



89, 



(0134) 
01 



81 



34 



^01^34 



(3.2) 



where Vj denotes the volume of the tetrahedron formed by all vertices except i, e.g. 4 — )• (0123), 



such that Vi = V( 



(0123)' 



Since 59, 



(0123) /r/)(0134) 



01 



) can only depend on Z23 (Z34 respectively), we can 



use equations (j3.2p in (|3.ip and obtain: 



81 



23 



^34 Vl 



81 IV- ^^-^^ 

0*134 i23 V2 

In general one finds (for a five vertex configuration with vanishing Cayley-Menger determinant) 
061117]: 

(3.4) 



8lij 




hkV-k 


dljk 







The actual sign depends on the geometric configuration under consideration. 

In addition to relation (j3.4p we need an analogous relation between deviations of edges not 
sharing a vertex. This can be derived from ()3.4p : Consider variations of three edge lengths lij, 
Ijk and Ikm such that wqi = 0. That is, lij can be understood as a function of Ijj. and Ikm- Then 



5lij 



Oljk + — OLkr, 



ULjk ^''km 

Now we restrict the variations further by requiring 5lij 
understand Ijk function of Ikm, that is Uj = I 







8lij 81 



'jk 



ij \ ljk{lkm)j Ik'i 

81 



n . (3.5) 

0, such that we have to additionally 
Thus one obtains for (13.51): 



81 jk 81 km 

8lij _ _ 

8lkm 



5lkm + 

lij Ijk 
Ijk Ikm 



8lkm 



5lkr, 



With ([331) we find (see also [17]) 



8lij 




IkmVk^m 


8lkm 




hjW] 



(3.6) 
(3.7) 

(3.8) 



To summarize (l3? 



3M: 



8lij 




IjkVk 




8lij 




Ikm^k^^ 


8ljk 






1 


8lkm 







where in fact p.4p is a special case of ()3.8p . 

In the following section we will use relations (|2.1ip . ([3 

9a, 
dl 



form %. 



and (|3.8p to compute terms of the 
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3.2 Computation of ^ 



The Hessian of the Regge action has entries of the form which we have to evaluate on 



dl ' 



configurations where uj^^^^^'^ = 0. As in the previous section we consider the initial configuration 
of the 3 — 2 move. We will start with the calculation of which is the derivative of the bulk 
deficit angle with respect to the bulk edge length. 



The equation of motion for the perturbations Aqi around the flat solution 1^^^ is given by 







E 



dlhdl 



01 



01 



■^01 



01 



01 



(3.9) 



Here h indicates edges in the boundary triangulation, and the sum is over all such edges. 

Now we know that these equation of motion specify Aqi , such that the linearized deficit angle 
at (01) is still flat. That is, if we choose the boundary perturbations, such that for instance 
only A23 ~ 5/23 is non- vanishing, we know that the ration of Aqi ~ and A23 ~ (5/23 has to 
satisfy (j3.8p . This specifies the ratio of the derivatives duoQi/dloi and dujQi/ dl2^. For the latter 
derivative, as ^23 is only included in one of the three tetrahedra we have via ()2.1ip 



(0123) 
01 



dl 



23 



^01^23 



This finally gives 



dujQi 




loiha 5^23 




'01 


VqVi 


dloi 




6V| 6I01 




6 





(3.10) 



(3.11) 



The actual sign is determined by the geometry and discussed in the next section. Note that we 
could have also used the lengths ?24 or Z34 instead of ^231 which would have however all lead to 
the same result. 

Next we consider terms of the form i.e. derivatives of the deficit angle at the bulk edge 
with respect to a boundary edge length. Note that for b = 23, 24, 34 the result is already given 
by (the analogue of) (|3.10|) . 

To find the derivative with respect to the remaining boundary lengths consider again ()3.9p 
with all boundary perturbations vanishing except, say Aoi ~ Sloi. Then, with the same line of 
arguments as used previously we can conclude 



01 



and hence 



dloi 




9^01 




dujQi 




loiJ-Oi 




dloi 




6 





(3.12) 



(3.13) 



where i £ {2,3,4} and j, k are such that i,j,k = 2,3,4. Again, the sign is determined by the 
geometry under consideration. 

Note that due to the symmetry of second derivatives of the Regge action we have 



duje 
dlpi 



Hence we can deduce terms of the form from . Thus only terms of the form remain 
to be computed, i.e. derivatives of exterior angles with respect to boundary edge lengths. 



(3.14) 
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To this end, remember that the initial configuration of the 3 — 2 Pachner move is flat. 
During the Pachner move the edge (01) is removed and replaced by a triangle (234), such that 
neither the intrinsic geometry (i. e. flatness) nor the extrinsic geometry (the embedding into flat 
space) of the boundary changes. In particular we will have that the extrinsic curvature angles 
io^ {h' ,loi{h')) = coincide in the initial and flnial configuration of the Pachner moves, 
involving three or two tetrahedra respectively. Here we understand Zqi a function of the 
boundary lengths If,/ as it is determined by the requirement of flatness. 

Now varying just one boundary edge lengths /;,', together with Zqi = loiih') as function of 
this lengths we obtain: 




(3.16) 



This gives flnally 



5a;f _ kh' WjVfnVn 



oh' oh' 6 Yln^n 

for b = (ij) and b' = (mn). Here s = ±1 denotes a sign, that will be determined in the next 
section. 



3.2.1 Determining the sign of 

In the previous section we have seen that in order to compute the full expression for the matrix 
elements of the Hessian, the actual sign of the derivatives of the form has to be determined. 

To be more precise, one only needs to determine the signs of where one has to treat the 
cases in which shares a vertex with Iqi and where it does not share a vertex separately. 

We start with the case where shares a vertex with Iqi, e.g. loi with i G {2,3,4}. In the 
derivation of the formula for ^j^, we considered variations of the edge lengths Iqi and loi, while 
keeping all other edge lengths fixed, under the condition that the triangulation is supposed to 
remain fiat. This allowed us to understand Iqi as a function of l^i, /qi = ^oi(^Oi)- To determine 
the sign of this dependence, consider Fig. [6) 

Assume that we enlarge hi slightly, i.e. 6hi > 0. If we do not change hi as well, the condition 
loqi = will be violated since all other edge lengths are fixed. However, if one allows hi to vary 
as well, the vertex (0) will be 'pushed' away from the vertex i, but since the edge lengths hj 
and hk are fixed, hi has to be increased, i.e. 6hi > 0. Hence: 

5/o. > =^ 6hi > ^ §^ = . (3.18) 

IMI 5hi hiVi 

We follow a similar line of argumentation for terms of the form with i,j £ {2,3,4}. 
Consider Fig. [Tj 

Assume that we slightly increase hj, i.e. 6hj > 0. Since the edge lenghts hk and Ijk are fixed, 
the vertex k is being 'pulled' towards the edge (ij). Furthermore the edge lengths hi, hj and 
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Figure 6: As one increases the edge lengths /qSj one also has to increase /qi in order to keep the 
triangulation flat, i.e. wqi = 





Figure 7: As one increases the edge lengths I23, one has to decrease Zqi in order to keep the 
triangulation flat, i.e. cjqi = 

/ofc are fixed, such that the vertex (0) is 'dragged' towards the edge (ij). This configuration can 
only remain flat, if Iqi is decreased, i.e. 6I01 < 0. Hence: 



5lij > 



6I01 < 



This is also consistent with (|3.1U|) . 



(3.19) 



3.3 Summary for the 3 — 2 move 

Let us summarize the results of the previous paragraphs. 

• In case either edge {ij) or the edge {km) is in the bulk, one obtains: 



H. 



(3) .__ 



du)ij 



dl 



km 



I K V-V-VLVrf, 



6 UnVn 



where 



1 ifie{0,l} 
else 



and the product in the denominator runs over all vertices in the triangulation. 



(3.20) 
(3.21) 
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• In case both the edges (ij) and (km) are in the boundary, one obtains: 

ivUkm)- Ql^^Ql^. Ql^^ y ) 6 rin^n Ohm ^' ' 

where denotes the exterior angle at the (boundary) edge {km) in the i tetrahedra 
configuration, Si is defined as above. 

Notice the simple form of the Hessian, 

^{ij),ikm) = Kij),(km) ^km) , (3-23) 

in particular that the second summand in ()3.23p factorizes. (Here H^.j^ ^^^^ = if either (ij) of 
(km) equals (01).) 

We have now all the prerequisites to discuss the (form-) invariance of the path integral 
associated to the 3 — 2 move. 

3.4 Invariance of the path integral 

For the 3 — 2 move we have to consider an expression of the following form: 



P3-2 = J dXoi fi{l) exp 



{ij),{km) 



(3.24) 



where: 



/x(Z) is a measure factor, which we assume to only depend on the background variables /, 
such that the configuration is flat. 

(3) 

{ij) {km) ^® (ij) (A;?7i)-matrix element of the Hessian in the three tetrahedra configura- 
tion, which we computed in the previous section, see ()3.20p . ()3.22p . 

Since Aqi is the only dynamical edge (variation) in the configuration under discussion, the 

^(3) 

'(01),(01) 

^'(01),(01) 



sign of ^(01) (01) crucial for the convergence of ()3.24p . In the sign convention introduces 
in equation ()2.2p H)^L > 0, such that ()3.24p converges. 



We can easily perform the integral in ()3.24p as it is a (partial) Gaussian integration. For an 
integral of the form 

/ = j dqi ... dqr exp Mgj (3.25) 

where M is a real, symmetric, positive-definite n x ?i-matrix and q = (qi) denotes a vector with 
i = 1, . . . ,r,r + 1, . . . n. Splitting the matrix M accordingly into submatrices 

M.(-»-J 

we can write 

/ = expl--ifU u\ with U := Uq - W^^ V . (3.27) 
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Using this result for the 3 — 2 move ()3.24p . we identify: 



H 



(3) 

{01),{01) 

{U{))b,b' 



dloi 



dly 



+ 



dlh dlh 



H, 



(2) 



doji 



'01 



Furthermore requiring form invariance of (j3.24p implies 

1 ^r(2) 



(3.28) 
(3.29) 

(3.30) 



P3_2 oc exp < 



(ij)7^(01),(fcm)^{01) 



2 ^{ij)ikm) ^km 



(3.31) 



such that in order to show that (|3.24p is invariant (on the level of the action), one has to show 
that U = H^'^\ which imphes 



dujQi dim 
dlb dly 



01 



dh 



01 



dlh dlh' 



(3.32) 



Note that we have aheady proven that ()3.32p holds due to the identity ()3.12p . This shows form 
invariance of the action. 

For the invariance of the measure ^{l) in ()3.24p we examine the contribution from the Gaus- 
sian integral: 

^/2^ 



^ ^ ' V (01), (01) 



dloi 



Vvi^ I V2V3V4 

VqVi 



01 



(3.33) 



Hence, choosing the measure factor as 



l,{l) 



Ue 



'VItx 



(3.34) 



we obtain a partition function, invariant under 3 — 2 Pachner moves. Here e denotes the edges 
and r the tetrahedra in the triangulation. 



3.5 4—1 move 

For 3D gravity, in addition to the 3 — 2 move, we have also to consider the 4 — 1 move. This 
move amounts to the subdivision of one tetrahedron, denoted by (1234), into four by adding 
one additional vertex (0), placing it inside the original tetrahedron and connecting it with all of 
the remaining vertices, see section P2.2.2I 

In contrast to the 3 — 2 move, the edge lengths of the new edges, i.e. the position of the new 
vertex inside the original tetrahedron, is not uniquely fixed. In fact the action is invariant under 
translations of the vertex (0) inside the tetrahedron (1234), such that one expects the Hessian 
matrix to have three null eigenvectors. In order to compute this matrix, terms of the form 
have to be evaluated just as in the 3 — 2 move. Following a similar derivation as in the previous 
section, one arrives at the following terms: 
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In case either the edge e = (ij) or edge e' = {km) are in the bulk, one obtains: 

_ ( iy,+s,+sk+s„,+i kjhm ^^j^kyfa ^^^^^^ 



(3.36) 



dlkmdhj 6 Hn 

where 

1 if i = 
else 

• In case both edges are in the boundary, one obtains: 

d'S ^ _ ^ (3.37) 

dkmdlij 6 rin^n (^^km 

where u^^^ denotes an exterior dihedral angle in the one tetrahedron configuration. 
Again, notice the simple form of the Hessian, 

= ^{ikikm) +^hij) hkm) , (3-38) 

with a factorizing summand. This form of the Hessian makes the appearance of null vectors 
obvious. 



3.5.1 Null eigenvectors 

'(Oi),(Oi) 



Since the pure bulk part H^^\ of the Hessian matrix factorizes, we can easily examine the 



condition for null vectors v: 

X] ^m,m^J = ^^Oi X] ^OjVj = (3.39) 
j j 

Hence, due to the factorizing form of the Hessian, we just have one condition for the null vectors. 
Therefore the Hessian has three null eigenvectors and of the four bulk degrees of freedom three 
are gauge. 

Furthermore we have to discuss the sign of the Hessian. The only non-vanishing eigen- 
value of the submatrix f/'(oi),{Oj) can also easily be determined due to the factorizing form to 
be X^j -ff(oj),(Oi)- This gives a negative eigenvalue, hence the Gaussian integral would not be 
convergent. This is a trace of the conformal mode problem in Euclidean gravity: the kinematic 
term of the conformal mode comes with the 'wrong' sign, so that the Euclidean action is not 
bounded from below. We see that the Pachner moves allow a nice isolation of this mode problem 
into the 4—1 moves. We will change the global sign for the action of the 4 — 1 move, such that 
the integral (over the one non-gauge mode) converges. This can be understood as selecting a 
complex contour for the integration for the conformal and the other modes separately, see [58] 
for a discussion in the continuum. 

3.5.2 Invariance of the path integral 

Similar to the 3 — 2 move, we have to consider 



4^1 



I n ^^0^ /^(/) exp I - Yl l^ltkikm) ^km \ (3.40) 



{ij),(km) 
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where ^{l) is again a measure factor, which we assume to depend only on background variables 
/, which have to make up a flat configuration. (j3.40p is again a partial Gaussian integral but 
with three gauge degrees of freedom, for which we will modify the general method of section [331 
Again the general form for the Gaussian integral is: 



I 



ji . . . dqr exp ^ —-(fMq 



(3.41) 



Since there are gauge degrees of freedom one integrates over the matrix M is singular. Assume 
that there are m gauge degrees of freedom such that we can split q in the following way: 



(gi, • • • 



1 Qr—mj Qr—m+l 5 ' 



=■■9 



, Qr,Qr+l, ■ ■ ■ , In) 



(3.42) 



(Here we assume that the transformation between g^-m+i, ■ ■ ■ ,qr and the m gauge parameters 
is not singular.) This implies the following split for the matrix M 



M 



Wo Vg V 

Vj Go Zo 
Zl Uo 



(3.43) 



where Wo is non-singular. Integrating out the degrees of freedom summarized in w one obtains: 



(r — m) 

(2vr)^ 
Vdet{Wo) 



exp |-i {fGg + fZu + ifZ^g + ifUu)^ (3.44) 



with 



G=Go-Vg'^Wo% , 
Z=Zo-VjW^^V , 
U =Uo - V^W^^V . 

Applying this formalism to the problem under discussion, one identifies (here i, j ^ 0,1): 



(3.45) 
(3.46) 
(3.47) 



(4) 

(01),{01) 



^b,b' + 



dio, 



(^o)(oi),{01) - ^ 

{Uo)b,b' 
iV)ioi),b = - 
We therefore obtain: 

(G')(oi),(oj) 

{Z){Oi),b 

{U)b,b' = 



doji 



01 



dloi 
dujoi dloi 
dlh dlhi 



'01 



duoi dloi 
dloi dloj 

duJoi dloi 
dloi dlh 



+ 



+ 



'01 



dloi 
dujoi 



(G'o)(oi),(oi) 

(^9)(01),(0i) 

iZo){oi),b 
diooi diooi 1IXT2I1 



dujoi dloi 
dloi dloj 
dujQi 

dloi 
duoi dloi 
dloi dlb 



dl 



01 



dloi dloj 
''^ diooi dujoi ilXT2t 







dloi dlh 







(1) 

b,b' 



dupi dloi 
dlh dlh' 



+ 



dtj, 



01 



dl 



01 



^ dujQi dioo 



'01 enj (1) 



dlh dlh' 



b,b' 



(3.48) 

(3.49) 
(3.50) 
(3.51) 
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This proves form invariance of the action, as the remaining term in the exponential corresponds 
to the action of the tetrahedron (1234) (after we have rotated back the global sign of the action). 

Note that after having only integrated over Aqi the other bulk variables Aoi, i = 2,3, 4 do 
not appear anymore in the exponential. 

Let us first consider how the measure factor is modified by the Gaussian integration over 
Aqi- The additional factor is given by 



Vdet(VFo) ^01 V ^0^1 

If we consider fj,{l) in ()3.40p to be the same measure which gives an invariant amplitude under 
the 3 — 2 move, namely 

n 

= ' (3.53) 

where e includes all boundary and bulk edges and r the four tetrahedra of the initial configura- 
tion, we obtain: 



4-1 ~ 



^ " 2^lh' h' > / y_ • (3.54) 

[ b,b' J ^ 



The remaining integral over the variables A02, • • • , A04 can be identified with an integration 
over the gauge orbit, which is given by the displacement of the inner vertex (0). As one can 
show [m HH] , see also appendix [Aj the following identity between integration measures holds 

d^x^ = '^y_ • (3.55) 

Here d^x^ is the integration measure of the Euclidean coordinates Xg, a = 1,2,3 of the vertex 
(0). The displacement of this vertex corresponds exactly to the gauge action of the discrete 
remnant of the diffeomorphisms [501 B5l [TO] . Hence we will replace the last factor in p.54p by 
1. This can be understood as resulting from a gauge fixing procedure, including the appropriate 
Faddeev- Popov determinant. (The numerical factors are chosen to conform with the integration 
measure found for the 3 — 2 move, however it is not possible to fix them uniquely.) 



4 Summary for 3D gravity 

For a general 3D triangulation we define the path integral for linearized Regge calculus by 

P:= /^^^ n dAeexp|-i^Fe,e'AeA,A . (4.1) 

le is the length of the edge e, Vr is the volume of the tetrahedron r, Ae is the edge length 
perturbation of the edge e and -ffg.e' is the e-e' matrix element of the Hessian matrix of the Regge 
action. The considerations conducted in the previous section show that (j4.ip is invariant under 
Pachner moves, in case one follows the gauge fixing and sign rotation procedure procedure for 
the 4 — 1 move discussed above. Hence (j4.ip does not depend on the choice of bulk triangulation 
and in this sense it is discretization independent. 
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Here we assigned the numerical pre-factor (127r)^^/^ to the edges of the triangulation (as 
the TT factors result from integration over edges). Another possibility would be to associate this 
pre-factor to the tetrahedra of the triangulation, in which case one needs to appropriately adjust 
the numerical constant in the gauge fixing prescription for the 4—1 move. 

Amazingly, the path integral measure which we found for linearized Regge calculus, coincides 
with the semi-classical limit of the Ponzano-Regge model \27\ \25\ [26] . This is a triangulation 
independenlH spin foam model for 3D quantum gravity. Here the numerical pre-factors (also 
given by (127r)~"'^/^) are associated to the tetrahedra. 

It would be interesting to see, whether this correspondence can be extended to 3D Regge 
calculus with a cosmological constant. This theory can be (classically) formulated in a triangu- 
lation independent way, by using curved tetrahedra [191 [39]. The corresponding quantization is 
given by the Turaev-Viro model [51], for which the semi-classical limit has been obtained |52j . 
Hence (|4.ip should give a triangulation independent amplitude for linearized Regge calculus 
with a (positive) cosmological constant by replacing Ig and Vr by their respective counterparts 
on the sphere, i.e. sin(/e), where le G [0,vr], and Vr, the determinant of the Gram matrix. 



5 Computation of the Hessian matrix in 4D 



We are now going to discuss the 4D case. We will proceed as for 3D, that is first determine the 
matrix elements of the Hessian and then consider the path integral for the Pachner moves. It 
will turn out that the 4 — 2 and 5 — 1 moves behave very similarly to the 3 — 2 and the 4 — 1 
move, respectively, in 3D. There is however an additional Pachner move in 4D, namely the 3 — 3, 
which is significantly different, and thus responsible for the non-trivial dynamics of 4D Regge 
gravity. 

The Hessian of the Regge action is given by 

d^S sr^ dAfituduJstAi, -siT^ d'^Afitu (bulk) d'^Astu (bdry) ,r-,\ 



dlijdlr, 



stu 



E 



stu sbulk 



dkidL 



-co 



^stu P)l..p)l "^stu 

stiiSbdry 

where ujstu is the deficit angle at the (bulk or boundary) triangle {stu). In the following we 
will not discuss the last two terms in (|5.ip because we will consider flat background solutions, 
i.e. (jj^^^'-^^ = 0. That is the second term in (jS.ip vanishes and the third term is unaffected by 
Pachner moves, since the extrinsic geometry, defined by the embedding into flat space, is not 
changed. (Furthermore this term is only multiplied by boundary perturbations, which are not 
integrated over in the path integrals.) 

Hence we define a reduced Hessian matrix: 



{ij),{km) 



A 



(5.2) 



which can be rewritten as the product of the following two matrices: 



H 



' dkn 

9^02 



dloi 



...\ 



/ 



dloi 
dloi 

V ; 



(5.3) 



Here the index summation is over all triangles in the triangulation. So as in 3D we have to 
compute terms of the form ^p-. To this end we will proceed similarly as in the 3D case, as 
described in the next sections. 



also requiring a gauge fixing procedure for tlie 4 — 1 moves 
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5.1 4 — 2 move 



5.1.1 Auxiliary formulas 

As in 3D we will need additional formulas to derive all entries of the Hessian matrix in a compact 
way. We start with the derivatives of the dihedral angles at a given triangle with respect to the 
length of the opposite edge (j2.1ip 



dir. 



12V 



where Aiji^ denotes the area of the triangle (ijk), Imn is the lengths of the edge (mn) and V is 
the volume of the 4-simplex (ijkmn). 

Now given a flat triangulation with six, i.e. D + 2, vertices, we will consider edge lengths 
variations of at least two edges under the condition that the triagulation remains flat, i.e. the 
deficit angles are vanishing, a; = 0. Then, along the same line of arguments as in 3D, we obtain 



In case the varied edges share a vertex: 



dlij 




Ijk Vk 


dljk 







(5.4) 



• In case the varied edges do not share a vertex: 



dlij 




I km Vj^Vm 


dlkm 







(5.5) 



where denotes the volume of the 4-simplex formed without the vertex k. To determine the 
sign of the derivatives in ()5.4l5.5p one has to consider the geometric set-up in detail. Note that 
(|5.4p and (jS.Sp are the exact 4D analogues of (j3.4p and p.Sp respectively. 

5.1.2 Computation of ^ 

Consider the 4 — 2 move, that is two 4-simplices (02345) and (12345), which share one common 
tetrahedron (2345). By connecting vertices (0) and (1), the two 4-simplices are split into four, 
namely (Olzj/c). This edge is the only bulk edge in the configuration with 4 simplices. As for the 
3 — 2 move we can use two facts, to specify the matrix elements of the Hessian. Namely on the 
one hand, that the equations of motions for the perturbation variable Aqi require flatness, on 
the other hand that lengths perturbations around flat space have to satisfy equations (j5.4|5.5p . 
The equation of motion is given by 

-f^(01),eAe + -H'(01)(01)^01 = . (5.6) 

e^(Ol) 

But as the perturbative solutions are also flat, the perturbation variables Ag have to satisfy the 
relations ()5.4|5.5p . Hence considering boundary data, such that only one Ag 7^ for e = [km) 7^ 
(01) and Aoi 7^ we can deduce 



H(Ql)^(km) 




5lm 




Ikm Vj^Vm 


^(01),(01) 




^^km 







(5.7) 
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To specify the (reduced) Hessian even further, we can use that the hnearized deficit angles 
also have to vanish. That is consider variations of two edge lengths, here J/qi = -^oi and 
(5^34 = A34. Using that the linearized deficit angle 610012 has to vanish we obtain 



6uJoi2 = 



dloi 

ot34 + — — 0/01 



012 



dl 



01 



^01^012 ^^34 



12^5 '^^01 



1531 



dloi 

^01^012 VgVj 



12 



(5.8) 
(5.9) 

(5.10) 
(5.11) 



The sign is again determined by the geometry under consideration. 

The rest of the derivations proceeds in the same way as for the 3 — 2 move. That is for the 
derivation of the bulk deficit angle with respect to a boundary edge length lb we obtain 



This can be repeated for all bulk deficit angles: 



duoii 




loiAoii y^y-i 


dloi 




12 W^Vn 



du}\ 



Oli 



dloi dlb 



(5.12) 



To determine the derivatives of the (boundary) exterior angles one again uses that these 
angles agree in both configurations of the 4 — 2 move. This results in 



(4) 
ijk 



dl 



and furthermore in 



doj, 



01 



(4) 
ijk 



12 



du. 



(2) 
ijk 



UnVn 



dloi 



(5.13) 



(5.14) 



dlb dlb dloi dlb 
The missing signs are dependent on the geometry under discussion and determined by similar 
considerations as in section 13.2.11 To summarize the results for the 4 — 2 move: 



In the case that either the triangle or the edge is in the bulk: 



^1 

81 



12 



where 



1 ifiG{0,l} 
else " 



(5.15) 



(5.16) 



In the case that both triangle and edge are in the boundary: 

ijk_ ^ ^ -^yi+Sj+Sk+s^+Sn+l ^mnAijk ^^jyiVniVn 



du^ 



dl 



12 



+ 



'^^ijk 

81 



(5.17) 



Note that as in 3D the formulas for ^ factorize. 
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5.1.3 Hessian matrix 

In order to complete the calculation for the (reduced) Hessian matrix, the terms ^f^^^ have 
to be summed up, where 

dAijk _ lij 2 . ,2 ;2 



f)].. 8 4., ihk + ljk-hj)- (5-18) 

- ^ij ; k 

Note that ()5.18p is only non-vanishing for four triangles in the triangulation for a given edge 
(ij). This implies (in case either (op) or (mn) are in the bulk): 

(4) _ ■\r^ dAjjk duijjk _ sr-^ dAppk duippk , . 

(ijk) k^o,p ^ 

" ^ 8 12 ^ ^ U,Vr 

symmetric in (op) <-> (mn) 



(5.20) 



where is a factor independent of the choice of (mn). In case (op) and (mn) are in the 
boundary, we obtain: 

(2) 

fT(4) _ r) I- 1 \So+Sr,+Sm+Sn ^op^mn VoVpVifiVn ST^ hp ^ ^^opfc 9^ 



(2) 

(op) , (m7i) 



ij(^) and 

are equal to matrices of second derivatives of the Regge action up to symmetric 
terms. Hence iif(^) and iif^^^ are also symmetric matrices. From this it follows that the factor 
Dop is the same for all choices of op: 

D:=Dop = Dmn . (5.22) 
Note that the elements H^^j^ ^^^^ satisfy the relation (j5.7p . 

5.1.4 Invariance of the path integral 

We have to consider: 



= j dXoi /x(0 exp I - l^itjUkm) ^km 



{ij),{km) 



(5.23) 



where the measure factor is supposed to depend only on the background edge lengths /. 
Note that (^^^ > such that (|5:23]l converges. (Note that in ()5.23p we did not include 

boundary terms which only depend on the boundary perturbations.) 

The computation is analogous to section [331 That is to show (form) invariance of the action, 
i. e. 



P4_2 oc exp < 



E lHf;mm)>^^J>^km\ (5.24) 
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we have to proof that 



h: 



(4) 

{ij),{mn) 



H, 



(4) 

(01),(01) 



E 

fc^O,l 



dloi dloi 



H 



(4) 



(4) 



(ij), (01)^^(01), (mn.) 



+ H, 



(2) 

{ij),{mn) 



dlij dloi dlmn 



E 



dAij^dujjjk 
dlij dl, 



01 



E 

> fc^O,l 



3/oi dlm^ 



(5.25) 



(ij) and (mn) denote two boundary edges. Applying ()5.14p to the last term on the right hand 
side of ()5.25p gives: 



E 



dloi dim. 



E 



dloi dl 



dl 



01 



01 



dlr 



which shows that (|5.25p holds. 

For the measure factor we examine the contribution from the Gaussian integral: 



(5.26) 



27r 



27r 



^/Wl 



V^V^V^V^ 1 



(4) 

(01),(01) 



ky^O,l dloi dloi 



'01 



VqVi 



D 



(5.27) 



Apart from the additional factor ()5.27p is of a similar form as ()3.33p for the 3 — 2 move in 
3D. Hence the invariant measure factor /x(/) should be proportional to: 



le 



V192^ 



UaWK 



(5.28) 



where e denotes the edges and A the 4-simplices in the triangulation. However, with this form, 
we will still get factors of by applying 4 — 2 Pachner moves. The factor D does not factorize 
into contributions that could be associated to 4-simplices or other subsimplices. It is rather a 
sum of terms involving the edge lengths of the entire triangulation associated to the 4 — 2 move. 
We will therefore defer the discussion of this factor D until after we have considered all the 
Pachner moves in 4D. 



5.2 5 — 1 move 

Let us now consider the 5 — 1 move in 4D. Again, many derivations will be similar to the ones for 
the 4 — 1 move in 3D. The 5 — 1 move corresponds to the subdivision of one 4-simplex, denoted 
by (12345), into five by adding one additional vertex (0), placing it inside the original 4-simplex 
and connecting it with all of the remaining vertices, see section 12.3.21 

Here, the edge lengths of the new edges, i.e. the position of the new vertex inside the original 
4-simplex, are not uniquely fixed. Accordingly there is a 4-parameter set of solutions and we 
expect to find four null modes in the Hessian. 

The derivation of the matrix elements for the Hessian proceeds as for the 4 — 2 move. We 
arrive at the following terms: 
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In case either the edge e = (op) or edge e' = {mn) are in the bulk, one obtains: 



:=D(5) 

where 

else ^^-^"^ 

and Fop-k is defined as in the previous section. Note, that as for the 4 — 2 move the factor 
D^^^ does not depend on the choice of indices (op) in ()5.29p . 

• In case both edges are in the boundary, one obtains: 

^(5) ^ y^+s,+s^+s„+l loplrnnVoVpVrnVn (^r,) (1) 

where H^^p^ („^„) denotes the {op)-{mn) matrix element of the Hessian of the one 4-simplex 
configuration. 

This gives all matrix elements of the Hessian of the 5 — 1 move. In the next section we will 
discuss the pure bulk terms more closely, in particular with respect to null eigenvectors. 

5.2.1 Null eigenvectors 

In this section we examine the pure bulk terms of the Hessian matrix, i.e. equation ()5.29p for 
edges (Oi) and (Oj) for arbitrary Then ()5.29p can be rewritten as: 

= _ (5) ^ ^ y ry ^(5) ^ (5 33) 

hoi hoj " ^ ' 



c 



(5) 

So, as in 3D, the bulk terms in the Hessian H^^.^ j^q^.^ factorize: 

^(S),(Oi) =c^OiS- • (5-33) 

Hence, following the argument in section [3.5.11 we can conclude that -f^(oi) (Oj) features four 
null vectors. The one non- vanishing eigenvalue is again given by i?(oj),(oj)) which amounts to 
a negative value. We will proceed as for the 4—1 move and change the global sign for the action 
associated to the 5 — 1 move. This can again be interpreted as taking care of the conformal 
factor problem in Euclidean gravity 



5.2.2 Invariance of the path integral 

Similar to the 4 — 2 move, we have to consider 



5-1 



/ n ^^0^ /^(O exp - l^llikn.) \ (5.34) 
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where ^{l) is a measure factor, which we assume to depend only on the background variables I, 
making up a flat configuration. (j5.34p is again a partial Gaussian integral but with four gauge 
degrees of freedom. The treatment of this integral will be completely analogous to the 4 — 1 
move. Due to the gauge modes we will first integrate only over Aqi- 

This integration will result in an exponential that is independent of the other variables 
Aoi, i = 2,3,4,5. To show that we obtain the action associated to the one remaining simplex 
(12345) we need to invoke the identities 



H. 



H 



(5) 

(01),(01) 



(5) 

{ij),{nin) 

E dAjjk dojjjk dim 



-1 



fr(5) fr(5) , o-(l) 

^fo-),(01)^(01),(mn) ^ ^iij),{mn) 



E 



dloi dl 
dim dl 



01 




(5.35) 



Similar to (|5.25p . we apply ()5.14p to the last term in ()5.35p . which gives: 



dim dlr. 



E 



dAmk duoik dl 



01 



dloi dim dl 



01 



(5.36) 



This proves (|5.35p and hence form invariance of the linearized action under the 5 — 1 move. 
The measure factor /_f(/) in (|5.34p changes by the Aqi integration by a factor 



Vl927r ViV^ViV^ 1 



V^(01),(01) ^ 



01 



^0^1 /D(5) 



(5.37) 



which turns out to be of a similar form as the contribution in the 4 — 2 move. If we consider 
fi{l) in ()5.34p to be given by 



I. 



we obtain for the path integral 

n 'fc 



5-1 



(5.38) 



b,b' 




hi 



^1 v^T92^ 



Vi 



(5.39) 



Also here the remaining integral can be identified with an integration over the gauge orbit, which 
is again given by the displacement of the inner vertex (0). In 4D we have the identity 
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where d'^x^ is the integration measure for the Euclidean coordinates Xq, a 
vertex (0). Hence we will replace the remaining integral in (|5.39p by 1. 



(5.40) 
1,2,3,4 of the 
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5.3 3 — 3 move 



In addition to the 5 — 1 and 4 — 2 moves, the set of Pachner moves in 4D includes the 3 — 3 
move. Here a complex of three 4-simplices is replaced with another complex of 3 4-simplices, 
such that the boundary triangulation is not changed. This move does not involve any bulk edge, 
hence, differing from all the other moves considered so far, we do not have an equation of motion 
associated to this move. 

There is another essential difference to the other Pachner moves, namely that the action is 
not invariant under 3 — 3 moves. Evaluating the (full) Regge action for the two configurations 
of the 3 — 3 move, one finds a difference, that grows quadratically with the deficit angle of the 
(only) bulk triangle |41j . Hence the action is not invariant in general under 3-3 moves, in fact, 
such an invariance applies only on flat configurations. This violation of the invariance of the 
action holds also for the quadratic action of the linearized theory, as can be expected from the 
behavior in the full theory and as can be checked explicitly on configurations with non-vanishing 
(linearized) curvature. 

The derivation of the (reduced) Hessian matrix for the 3 — 3 move proceeds in a slightly 
different way, as we now have to take into account that the boundary perturbations might 
describe curvature. The result will however have the same structure as for the other Pachner 
moves. 

To start the derivation note that in both configurations A (with simplices (01234), (01235), 
(01245)) and B (with simplices (01345), (02345), (12345) there is only one bulk triangle, namely 
(012) and (345) respectively in A and B. The vanishing of the linearized deficit angles defines 
boundary perturbations in flat directions: 

(ij) (ij) 

For such flat variations Xij we can again derive, along the same arguments as in section [5. l-H 
the relations (15.51) 

(5.42) 









Afcm 







Note that equation ()5.4ip also implies that the gradients of the deficit angles in the two con- 
figurations are parallel to each other. That is the space of flat boundary perturbations Xij in 
both conflgurations is the same, the linearized curvature will however have different values in 
the general case, if evaluated on the same set of (non-flat) boundary perturbations Xij. 
Now, starting with the derivatives 

9^012 ^ ^^012^^ ^ ^45^123 ^^^45 ^ ^^345^^ ^ ^12^345 A'i) 

dki dlu 12^5 ' dh2 dh2 12^0 ' 

and using (|5.41|) and (|5.42|) . assuming that only two appropriately chosen length perturbations 
Xij,Xmn are not vanishing, we can obtain all other derivatives of the bulk deficit angle in both 
configurations. The signs can again be determined as in section Hi. 2. II 
The result is given by 

'9w^i2 kjAou y-^y\^-2W-j 



dli, ' ' 12 W^V^ 

dk, ^ 12 w.y, ^'-^'^ 
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Here we defined the sign factors as 



sf-\' ■"^*°-^!' . sf^il -f'^P-^l . (5.5) 



else ' * [0 else 

With the understanding that uj^^r^ = = ujqi2 '^^^ write the relation between the deriva- 
tives as 

^^opk _ ^^^^,_.yf+sf+st+s^+s^ kjAppk VoVpVkViV-j 

dkj dk, l—L ^ ' 12 Dp^P ^ ■ ^ 

where {opk) is either the set (012) or (345). Note that (|5.46p is consistent with ()5.45p under 
change A -H- i.e. change of sign. 

We will soon discover that (j5.46p holds also for the other (boundary) angles. To this end we 
use that the linearized boundary extrinsic curvature angles coincide in both configurations A 
and if evaluated on flat boundary perturbations. Hence we can conclude that the difference 
of the gradients of a given boundary angle has to be proportional to the gradient of one of the 
bulk angles, i. e. 



Again we can start with an especially simple derivative, i.e. 

(5.48) 



_ h2A345 V^ViV^ViV^ dug^ 



dlu 12 UpVp ' dl 



12 



to get hold of all the other derivatives of this exterior curvature angle. In this way we obtain 

for all the boundary and bulk angles. 

To finally arrive at the (reduced) Hessian, we have to multiply this result ()5.49p with the 
area derivatives as in (j5.19p . This allows us to express the difference of the (reduced) Hessians 
in the A and B configurations as 

H^opUmn)-H(^op),(rnn) " [-^) ]-[^ y_ 2^ (^_lr^ ^op;fc^fc (5.50) 



Also here D does not depend on the choice of indices o,p in (1^301) . Note that = -D^ . 
Unless D"^ = = 0, the (quadratic) action of the linearized theory is different in the A and 
the B configuration. This equality D"^ = = does not hold on general (flat) background 
configurations, but might hold in very symmetric cases. 

Furthermore, a measure of the form 



Kl) = ^ (5.51) 

is only invariant under the 3 — 3 move in the case that 

VoViV2 = V^ViV^ . (5.52) 
Again, this equality does not hold for generic cases. 
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6 Summary for 4D gravity 



4D classical Regge calculus is invariant under the 4 — 2 and 5 — 1 moves, but not under the 3 — 3 
moves. Our calculations provided the evidence for the linearized theory, in particular isolating 
the invariance breaking term for the 3 — 3 move. 

But this invariance behavior also holds for the full theory: there is always a flat solution to 
the equation of motions associated to the 4 — 2 and 5 — 1 Pachner moves. Hence the contribution 
from the bulk to the Hamilton-Jacobi function is vanishing. The Hamilton-Jacobi function is 
therefore just given by the boundary terms, which do not change under the 4 — 2 and 5 — 1 
moves. 

Concerning the quantum theory for linearized Regge Calculus in 4D, we define the path 
integral for general triangulations as 

P:= rr dAeexp|-ii?,,e'AeAe4 • (6.1) 

le denotes the length of edge e, Va the volume of 4-simplex A, Ag is the edge length perturbation 
of edge e and -ffg.e' is the e-e' matrix element of the Hessian matrix of the Regge action. 

In the previous section we have shown that ()6.ip is invariant ~ modulo the factor D ~ under 

4 — 2 and 5 — 1 Pachner moves (using the gauge fixing conventions discussed above), but in 
general not under the 3 — 3 Pachner move. The non-invariance under 3 — 3 moves is already 
present in the classical theory and should be overcome by constructing a perfect discretization 

It might be possible to implement a full invariance of the path integral under either the 4 — 2 
or the 5 — 1 move, that is by including the factor \/TD into the measure. For the 4 — 2 moves 
one would need to associate a corresponding factor to the edges of the triangulation, for the 

5 — 1 move rather to the vertices. (Alternatively, one would have to change the gauge fixing 
procedure for the 5 — 1 move, i.e. the factor associated to the gauge orbit, but this seems to be 
rather unnatural.) Still there are several open questions left to address, as how to generalize the 
definition of the D factors to more complicated triangulations ( the bulk edges in the Pachner 
moves are always shared by four triangles) and how the D factors associated to boundary edges 
or vertices will interfere. Furthermore, the factor D is slightly non-local, but its actual form 
might be due to the linearized theory. 

Here it might be helpful to reconsider the topological BF theory, from which gravity can be 
obtained by implementing (simplicity) constraints. This is the route followed by spin foams. The 
advantage of this approach is, that a triangulation invariant path integral can be constructed 
for BF theory. To apply this to Regge calculus one would need a formulation based on the 
same geometric variables as used in 4D BF theory. Such a formulation is provided by area- 
angle Regge calculus [38]. The corresponding action can also be split into a piece describing a 
topological theory and constraints acting in the same way as the simplicity constraints. Studying 
this action might help to construct a triangulation independent quantum theory describing flat 
space dynamics. For other work in this direction, related to BF theory see [54: \ 155 1 [56]. 

Despite all these subtleties and drawbacks, the simple form of the Hessian matrix for all 
Pachner moves and its similar form to the 3D case are remarkable. Therefore it will be very 
interesting to compare our results to spin foam asymptotics and possibly help to fix measure 
ambiguities (by requiring invariance under Pachner moves) there. 
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7 Discussion 



In this work we provided extensive analytical calculations for linearized Regge calculus, very 
much enlightening the structure of the theory. In particular we obtained the linearized Regge 
actions associated to all the Pachner moves in 3D and 4D, explicitly showing that the Regge 
actiorlfl is invariant under all Pachner moves, with the exception of the 3 — 3 move. We isolated 
the gauge symmetries and the conformal factor problem, which both are potential sources for 
divergencies. Amazingly the structure of the linearized Regge actions associated to the Pachner 
moves lead in all cases to a very transparent factorizing structure, similar in 3D and in 4D. 
These formulae might be also helpful in other contexts, for instance in a canonical formulation 
of Regge calculus [52] or in numerical larger scale calculations. 

Furthermore we proposed a dynamical principle to fix the measure for Regge calculus, namely 
to consider the behavior of the theory under Pachner moves. Restricting to a local ansatz as 
in ()1.2p this fixes the measure uniquely. Indeed the invariance under change of triangulation is 
related to an implementation of diffeomorphism symmetry \W\ \T9\ \TE[ [T5] . This condition can 
therefore be understood as requiring an anomaly-free measure, which can be expected to be 
unique. A simple reason for this is that for a theory completely invariant under changes of the 
triangulation, there is also no (bulk) discretization scale. That is the only discretization scale is 
provided by the boundary data. For compact manifolds without boundary the continuum limit 
is even trivial, as such a limit would be obtained via a refinement of the triangulation [T7j- In 
other words a triangulation independent path integral provides already the continuum result. 

This is the reason why we cannot expect to obtain a fully triangulation independent local 
theory in 4D. However one might ask for invariance of the quantum theory under the same set 
of local triangulation changes under which the local classical theory is invariant. Such a set can 
be understood as trivial subdivisions of the triangulation, as the associated equations of motions 
lead to flat space-time. The question therefore is, whether one can define a topological sub- 
sector of the theory [6l[55], which would provide a quantum description of flat space dynamics, 
see also [56j- Such an invariance of the theory under trivial subdivisions seems also be crucial to 
realize scenarios as proposed in |17] . i. e. the convergence of the theory to a topological sector, 
under refinement. 

We found a path integral measure for linearized Regge calculus, which provides such an 
invariance, in 4D modulo a factor, which features a certain non-local structure. In 3D we found 
an exactly invariant measure, which also coincides with the asymptotics of the (triangulation 
independent) Ponzano Regge model. This is quite astonishing, as we performed a calculation 
in the linearized theory. Furthermore the Ponzano Regge model includes in addition also a sum 
over orientations, which we do not consider here. The question arises, whether this result can 
be extended to the full non-linear theory and shed light on the problem, whether to include a 
sum over orientations into a quantum gravity path integral or not \57\ I58|. 

The factor appearing in 4D, disturbing invariance at least under the 5 — 1 and 4 — 2 moves, 
is related to a transformation from area to length variables. It might therefore be helpful, 
in order to further enlighten this issue, to consider area-angle Regge calculus |38[ 139] . This 
formulation allows a split into a topological theory, which would be triangulation independent, 
and constraints. Another possibility to obtain path integral measures is a derivation from the 
canonical theory, which has recently became available for Regge calculus [l2]. Indeed the path 
integral measure is important to obtain correlation functions, which are annihilated by the 
Hamiltonian constraints [16] . 

''This invariance result holds also for the full theory. 
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A fully triangulation independent theory can be constructed via the method of perfect dis- 
cretizations [15^ [T9l [THl I16j . which is based on a Wilsonian renormalization flow. This has the 
advantage of providing at the same time informations on the continuum limit of the theory. 
Here two different strategies can be thought of. One is based on local considerations, namely 
to study the behavior of a given theory under local refinements, e. g. Pachner moves, see for 
instance [31] for related studies in (topological) spin foam models. This can result in recursion 
relations, whose fixed points provide the continuum limit (and perfect discretization) of the path 
integral, see [16] for an example in ID. Another strategy is to extract the large scale behavior, 
which might depend on the choice of measure [131 . First steps towards extracting large scale 
behavior of (simplified) spin foam models via real space renormalization can be found in |60j . 
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A Euclidean integration measure 

The usual Lebesgue measure of D-dimensional Euclidean space can be rewritten with respect 
to the edge lengths of a (non-degenerate) L)-simplex |49j . 

Assume D + 1 vertices making up a D-simplex embedded in M^, their positions given by 
{xi}j=o,...,D) such that the /^-simplex is not degenerate, i.e. its D-volume is non-vanishing. 
Next, we define the position of the vertices of the D-simplex with respect to one of its vertices 
by defining li := x-i — xq. Since the D-simplex is non-degenerate, the set of vectors {/i}i=i,...,D 
form a (non-orthonormal) basis of M^, where the lengths of the vectors li give the edge lengths 
loi of the D-simplex. To write the Lebesgue measure in these coordinates, one has to compute 
the Jacobian of the linear function which maps the orthonormal basis {ei} to {li}. To simplify 
notation, we will denote y := xq. 

D D /,9x D . ,^ ^ .o\ D 



W loi dloi 



i=l ^ 



[y-XiY 



Wdyi 



det 



d {y 



Xi 



dyj 



(A.l) 



=1 i=l ^ ' i=l ^ ' i=\ 

where the determinant in the last term is the Jacobian of the coordinate transformation. For 
the matrix elements of the Jacobian one obtains: 

^2 



d {y 



-3 



■ {y- Xi) 



dy^ 2 

Given (|A.2p . the Jacobian can be rewritten in terms of the volume of the D-simplex: 



|det (ej • {y-xi))\ = ^det {{y-Xj) ■ {y - Xi)) = D\V . 
Using (|A.2IA.3|) in (IXH) gives: 



D 



D 



J{lQ.dlQi = J{dyiD\V 



i=l 



D 



i=l 



Y{dyi 



i=l 



Ui=i ki d l 
D\V 



Oi 



(A.2) 
(A.3) 

(A.4) 
(A.5) 
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Prom (|A.5P one obtains ([335]) in 3D and ([5^ in 4D. 



B Determinant formula for D 



In the calculations of the Hessian matrix in 4D, see section O we have encountered the factors Dop 
(see for instance ()5.19p ) which are slightly non-local. In this section we will present a different 
way to compute these factors in terms of a determinant of a matrix which will additionally 
allow us to show that (some) factors Dop are equal without using that the Hessian matrix is 
symmetric. 

Consider six vertices embedded in making up a triangulation which can be modified by 
one of the Pachner moves discussed in section [231 We define the position vectors of the vertices 
with respect to vertex (0), which we place in the origin of the coordinate system for simplicity. 
Hence the position vector of vertex (i) is defined by x* = x* — x*', its components are denoted 
by x\ with a = 1,2,3,4. Given this definition, consider the following determinant: 



det 



(fl)2 (f2)2 



(x5)2 



(B.l) 



We will show that this determinant is proportional to the factors Doii for z = 1, . . . , 5. 

The determinant of a matrix remains unchanged if a scalar multiple of one of its rows / 
columns is added to one of its rows / columns respectively. Hence we will subtract x\ times the 
ath row from the last for 1 < a < 5. Then one obtains for instance for i = 1: 



det 



X2 
2 



Xa 



->2 —1 
X ■ X 



(f2)2 

The terms in the last row can be rewritten as 



'5^2 ^5 ^ ^1 



{x^y -x''-x' = -{{xy - {x'Y + (x" 



-/2 _/2 1/2 
''■Ok 'Ol^'lfe 



(B.2) 



(B.3) 



where -Foi;fc is the same factor as in (|5.18p . Expanding ()B.2p with respect to the last row, 
the subdeterminants of the matrix correspond to (oriented) volumes of 4-simplices (see (|A.3p ). 
Hence one obtains: 



det[...] = 4! J^(-l)^-isiFoi;,y. oc D 



01 



(B.4) 



i=2 



where Si are appropriate sign factors taking into account the orientation of the volumes. By 
analogous considerations, one computes the other factors D^j, j = 2,3,4,5, from the same 
determinant, which implies that these factors are all equal (although the symmetry of the Hessian 
matrix has not been used). The signs Si depend on the relative orientation of the vectors Xj of 
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the respective 4-simplex and hence they depend on the geometry under discussion. We would 
hke to demonstrate that exphcitly with the example of the 4 — 2 move. 

Consider the Dqi for the 4 — 2 move, see ()5.19p . Four of the five edge vectors Xi of the 
vertices, defined with respect to the vertex (0), form a basis of since the 4-simplices are not 
degenerated. In the following we assume that X2,---,X5 is a positively oriented orthonormal 
basis. Hence, we can rewrite the determinant (|B.2p as 

YA=2'"i^i ^2 X2, Xi X5 
Foi-2 Foi;3 Foi-A Foi;5 



det 



(B.5) 



where the Vi are the coefficients of xi in the basis formed by X2, • • • , a?5- Note that the configu- 
ration of the 4 — 2 move can be chosen such that Vi > 0. If (|B.5|) is expanded with respect to 
the last row, the submatrices are of the following form: 

5 



det 



^ijkl "^i ) 



(B.6) 



^ ^ fjXj Xj Xfc Xl 
1=2 

where eij^i is the Levi-Cevita symbol in 4D. Hence the signs Sj in (|B.4p are alternating, which, 
together with the alternating signs (—1)*"^ verifies formula (|B.4p . 
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